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ABSTRACT 


A  computer  program  is  described  which  enables  on-line  display  of  chosen 
sections  of  time-series  data  files.  Files  can  be  stored  either  on  magnetic 
tape  or  on  disk,  and  two  files  from  the  same  unit  can  be  displayed 
simultaneously.  The  program  is  fully  interactive,  flexible,  and  simple  to 
use.  Default  options  are  used  throughout  whenever  possible.  Fourier 
transforms  can  also  be  computed,  and  the  resultant  power  spectra  plotted. 
Examples  are  presented  of  both  the  output  display  and  program  usage. 
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1 .  INTRODUCTION 


la  experiments  where  data  are  recorded  either  continuously  or  for  long 
periods  of  time,  it  is  fairly  clear  that  vast  amounts  of  information  can  soon 

be  accumulated.  Only  in  a  few  cases,  such  as  perhaps  the  study  of  some 

ambient  property,  is  analysis  of  all  the  available  data  required.  In  most 

cases,  only  small  sections  of  the  data  are  relevant,  but  the  remainder  must 

be  recorded  to  ensure  that  the  significant  events  are  themselves  logged. 
This  is  particularly  true  when  studying  relatively  infrequent  short-lived 
natural  phenomena  whose  times  of  occurrence  are  not  known  in  advance. 

There  is  thus  a  need  for  a  means  of  readily  and  speedily  isolating  the 
comparatively  short  data  sections  containing  the  'relevant  events'  and 
discriminating  against  the  much  longer  but  less  significant  background  data. 

Batch  processing  is  not  well  suited  for  such  selection,  since  any  one 
particular  portion  of  data  may  have  to  be  viewed  under  several  different 
conditions  to  decide  if  it  is  really  of  use.  Furthermore,  in  more  basic 
research,  there  may  well  be  little  preconceived  idea  of  the  form  of  the  event 
for  which  one  is  searching,  thus  making  the  development  of  automatic  scanning 
algorithms  impracticable.  On-line  operation,  on  the  other  hand,  is  ideal 
for  such  previewing  and  hence  for  the  selection  of  data  segments  that  merit 
more  detailed  analysis.  An  additional  benefit,  apart  from  the  obvious 
saving  in  time,  is  that  the  immediate  response  gives  the  user  a  much  better 
feel  for  his  data,  and  thus  aids  in  his  selection. 

The  purpose  of  the  present  article  therefore,  is  to  present  a  computer 
program  that  enables  such  on-line  previewing,  and  which  we  have  found  useful 
in  our  own  investigations  of  long  term  geophysically  related  data.  The  user 
can  move  at  will  to  different  sections  of  his  data,  replot  interesting 
sections  using  different  scale  conditions,  and  as  a  further  aid,  request  a 
Fourier  transform  and  view  the  resultant  power  spectrum.  To  aid  in 

comparison  between  different  measurements,  two  time  series  can  be  manipulated 
simultaneously.  The  program  is  currently  being  run  using  a  Tektronix  CRT 
terminal  connected  to  a  Sigma  7  computer  operating  under  CP-V.  Little 


difficulty  is  foreseen  in  transferring  the  program  to  an  alternative  system 
that  supports  Calcoraptype  calls  for  on-line  video  terminals.  Areas  that  are 
system  dependent  have  been  placed  in  subroutines  and  are  discussed  below. 

The  following  sections  discuss  the  use  of  the  program  and  describe  the 
routines  of  TEKPLOT.  ttt  the  appendix  we  give  flow-charts  for  the  two 
largest  routines  and  give  a  listing  of  the  complete  program. 

2 .  PROGRAM  USAGE 

Before  discussing  the  program  itself,  it  is  convenient  to  give  an 
overview  of  the  use  of  TEKPLOT.  As  will  be  seen,  the  program  has  been 
designed  for  a  minimum  of  user  input,  and  consequently,  a  number  of  default 
values  have  been  set  up.  Provided  the  user  is  satisfied  with  the  individual 
defaults,  he  need  only  enter  a  carriage  return  (hereafter  CR)  in  response  to 
several  questions.  Certain  other  options  can  be  ignored,  and  scaling 
conditions  chosen  automatically,  by  simple  CR  responses.  Sample  input  and 
plots  are  shown  in  Figures  1-3,  and  a  list  of  the  necessary  input  parameters, 
together  with  their  default  values,  is  given  in  Table  1. 

2.1.  INITIALISATION 

2.1.1.  FILE  SPECIFICATION 

After  starting  the  program,  the  user  is  prompted  for  the  name  of  the 
file  from  which  data  are  to  be  read  and  for  the  number  and  label  of  the  tape 
on  which  the  file  resides.  Entering  9998  or  9999  at  this  point*  informs 
the  program  that  the  file  resides  on  disk.  CR  defaults  to  disk  (9998). 

The  user  is  then  prompted  for  a  second  file  which  is  assumed  to  be  on 
the  same  unit.  If  only  one  file  is  desired,  then  CR  is  entered. 

t  In  the  present  system,  the  number  9998  refers  to  a  disk  used  for  permanent 
storage,  whereas  9999  indicates  a  temporary  storage  disk 


FILE 


X-01H 


TAPE  9996 


itUEF  FILE  NAME  IN  2A4 
TX-01N 

ENTER  Tape  .  14,  3A4,  CR  DEFAULTS  TO  9998 

FILE  2  CR  *  NONE 
’XTESTT 

ENTER  NSTART,  TINE  LENGTH,  ANO  SECS  PER  PAGE 
DEFAULTS  ARE  1,  28-  AND  10  FORHAT  110 
•30. 

’12. 

’12. 

ENTEP  SAHPLE  PATE  IN  15.  CR  IS  40  SANP  PER  SEC 
’45- 

ENTER  STAPT  TINE  OF  TAPE  CR  TO  IGNORE  THIS  OPTION  12.213 
’12  30  05 

ENTER  SCALING  CR  FOR  PER  PACE,  t  FOR  NMOLE  DISPLAY,  2  TO  BE  INPUT 

SCALE  PER  PAGE 
TAPE  9998  FILE  X-01H 

XTESTT 

TAPE  START  12  30  5 

START  AT  30  NUMBER  OF  POINTS  540 

SECONDS  PER  PAGE  12  DATA  POINTS  PER  INCH  54  00 


ENTER  CR  TO  CONTINUE.  A  TO  ABORT  ?■ 


Figure  2 


An  Example  of  the  Time  Series  Display 
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Figure  3.  An  Example  of  the  Power  Spectra  Display. 
Table  1.  Input  Parameters  and  Default  Values 


Variable 


Purpose 


Default 


IFILE 

name  of  first  file 

NTAPE 

tape  number  and  label 

IFILE2 

second  file  name 

NSTART 

starting  point  in  file 

NDISP 

total  number  of  seconds 
to  display 

NSEC 

seconds  to  display  per 
page 

ISP 

sample  rate  per  second 

IHST,  IMST ,  ISST 

start  time  of  tape 

ISCALE 

scaling  mode 

9998  (disk) 
option  ignored 
1 

20 

10 

40 

option  ignored 
0 


additional  parameters  whenever  used 

YS  range  per  division  for 

time  series  plots 

PLTMIN  minimum  for  power  spectra 

plot 

range  per  division  (in  dB) 
for  plot  of  spectra  plot 


PLTSTP 
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2.1.2.  DATA  SEGMENT  SPECIFICATION 

The  next  four  parameters  refer  respectively  to  the  point  in  the  file 
from  which  to  start  plotting,  the  total  number  of  seconds  of  data  to  be 
plotted,  the  number  of  seconds  of  data  to  be  plotted  per  page,  and  the 
sampling  rate  of  the  data  in  points  per  second  (which  must  be  integral).  CR 
response  to  any  of  these  prompts  results  in  the  displayed  default  value  being 
subst ituted . 

2.1.3.  TIME  OPTION 

To  the  next  prompt,  concerning  the  time  option,  the  starting  time  of 
the  data  file  can  be  entered.  Thereafter,  each  page  of  the  plot  will 
indicate  the  time  of  the  start  of  the  particular  segment  being  displayed.  If 
time  keeping  is  not  required,  simply  enter  CR. 

2.1. A.  TIME  SERIES  SCALING  OPTIONS 

In  response  to  the  next  prompt,  three  alternatives  can  be  given  which 
affect  the  manner  in  which  the  time  series  is  scaled  before  plotting. 

In  mode  0  scaling  (CR  entered)  each  page  is  treated  independently,  and 
the  scaling  conditions  are  chosen  automatically  so  that  the  plot  remains  in 
the  range  of  the  screen.  Furthermore,  if  two  files  are  being  plotted,  each 
time  series  is  scaled  independently.  This  mode  is  the  most  useful  when  the 
two  time  series  have  markedly  different  amplitudes,  and  when  the  data  have  a 
wide  dynamic  range. 

In  mode  1  scaling  (a  response  of  1)  the  maximum  and  minimum  values  are 
found  for  the  whole  of  the  file  one  data  segment  defined  in  2.1.2.  These 
values  are  then  used  to  determine  the  upper  and  lower  edges  of  every  time 
series  plot.  The  second  file  is  plotted  using  the  same  scaling  factors. 
This  mode  is  used  when  many  pages  of  a  time  series  have  to  be  used  with  the 
same  scaling,  or  when  a  quick  visual  comparison  of  two  files  is  required. 
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In  mode  2  scaling  the  user  is  prompted  to  enter  the  value  appropriate 
to  the  ordinate  division.  The  program  then  positions  the  plot  automatically 
so  that  the  minimum  of  the  plot  lies  at  the  bottom  of  the  y-axis.  This 
positioning  is  done  independently  for  each  page,  and  for  each  plot  on  the 
page  if  two  files  are  being  plotted.  (  It  should  be  noted  that  when  two 
time  series  are  being  plotted,  then  the  scale  factor  input  is  doubled  before 
use. ) 


2.2.  POST  PLOT  INPUT 

After  the  above  parameters  are  entered,  the  user  is  given  a  pause  to 
review  his  input.  Typing  CR  instigates  the  plot,  whereas  entering  'A' 

returns  the  program  to  2.1.2  and  requests  a  new  time  series. 

When  each  page  of  the  time  series  plot  is  complete,  the  user  has  four 
choices.  Entering  CR  causes  the  program  to  continue  with  the  next  page  of 
the  time  series,  whereas  the  'A'  option  causes  all  remaining  plots  of  the 
time  series  to  be  ignored  and  the  program  to  return  to  2.1.2.  The  'F1 
option  causes  a  Fourier  Transform  to  be  computed,  and  the  power  spectrum 
plotted  using  default  scaling  factors.  Power  is  plotted  in  dB,  and  the 
frequency  range  is  from  d.c.  to  one  half  of  the  sampling  frequency. 

The  'S'  option  also  plots  the  power  spectrum,  but  the  user  is  first 
prompted  for  the  minimum  and  the  range  per  division  of  the  power  axis. 
These  values  then  become  the  default  values  for  subsequent  * F *  calls  and 
remain  unchanged  until  a  further  'S'  call. 

2.3.  POST  POWER  SPECTRA  INPUT 

On  completion  of  a  power  spectrum  plot  the  program  pauses.  Entry  of 
CR  at  this  point  causes  the  program  to  move  to  the  next  page  of  the  time 
series.  Entry  of  'A'  returns  the  program  to  2.1.2. 

When  all  pages  of  the  specified  part  of  the  time  series  are  exhausted, 
the  program  prompts  for  a  further  segment  (section  2.1.2).  If  no  further 
plots  are  required,  and  end  of  file  (ESC  F)  will  terminate  the  run. 


3.  THE  FOURIER  TRANSFORM 


The  Fourier  Transform  routine  invoked  is  a  conventional  single  sided 
fast  transform.  Before  transforming,  the  mean  is  subtracted  from  the  data, 
and  the  data  tapered  by  a  sine  squared  function.  This  tapering  is  applied 
over  a  range  of  10%  of  the  data,  taken  at  each  end  of  the  data  (see  Figure 
4).  Zeroes  are  appended  whenever  necessary  to  increase  the  number  of  data 
points  to  an  integral  power  of  two. 


10%  OF  NDO 


Figure  4.  The  tapering  applied  in  the  POWSPC  routine,  a)  shows  the  sine 
squared  window,  and  b)  shows  a  cosine  data  train,  of  length  240 
elements,  after  tapering. 


4.  PROGRAM  DESCRIPTION 


In  this  section  we  present  a  description  of  the  various  routines  used 
in  TFKPLOT.  Supplementary  description  can  be  found  in  the  flow  charts  and 
listing  of  Appendix  A.  Any  areas  of  subroutines  that  may  contain  system 
dependent  code  are  identified. 

4.1.  MAIN  PROGRAM 

Two  of  the  principal  arrays,  DATA  and  DAT2,  are  used  to  contain  data 
appropriate  to  one  page  of  the  time  series  for  file  one  and  file  two 
respectively.  The  array  Xl  is  used  to  contain  the  generated  abscissa  values 
needed  for  plotting.  These  arrays  are  dimensioned  NMAX  +  2,  where  NMAX  is 
the  largest  number  of  points  per  page  that  the  program  is  configured  to  use. 
For  our  application  an  NMAX  of  4096  was  satisfactory,  but  this  could  be 
increased  provided  sufficient  core  is  available. 

After  opening  the  plotting  routines  and  clearing  the  screen,  the 
program  prompts  the  user  for  pertinent  information  concerning  the  time  series 
to  be  plotted.  The  responses  are  subsequently  analysed  for  CU  default 
replies,  and  the  necessary  default  parameters  set  up. 

Various  other  parameters  and  control  variables  are  then  set  up.  These 
include  three  logical  variables,  DOTIME,  TWOPLT  and  DISK,  which  indicate  that 
the  time  option  is  being  used,  that  two  files  are  being  analysed 
simultaneously,  and  that  the  original  data  files  were  stored  on  disk. 

Parameter  set  up  is  followed  bv  relevant  parameter  printout,  and  then  a 
pause  to  enable  the  user  to  check  the  variables.  A  CR  response  causes  the 
program  to  continue,  whereas  an  'A'  response  causes  the  program  to  return  to 
statement  number  884  and  request  the  input  parameters  once  again. 


9 


The  program  then  checks  to  see  if  two  files  are  being  used,  and  if  so, 
whether  they  reside  on  disk.  If  the  files  are  found  to  be  on  magnetic  tape, 
then  file  two  is  first  copied  to  disk.  This  copy  prevents  multiple 
rewinding  of  the  tape  during  the  plotting  sequence. 

The  program  then  checks  the  scale  mode,  and  if  mode  one  is  being  used, 
then  the  maximum  and  minimum  of  the  first  time  series  are  found,  and  the 
appropriate  scaling  factors  determined  for  later  use.  During  this  process, 
if  file  one  was  stored  on  tape,  the  time  series  is  copied  to  disk.  This 
markedly  increases  efficiency  by  avoiding  duplicate  tape  input  during  the 
plotting  loop. 

The  program  then  enters  the  main  plotting  loop.  Provided  that  the 
previous  page  was  not  the  last  page  of  the  specified  time  series,  a  new  page 
of  data  is  read  into  the  array  DATA,  and  if  necessary  into  DAT2.  The 
abscissa  values  are  computed  and  stored  in  the  array  XI.  The  scale  mode  is 
again  checked,  and  the  minimum  of  the  plot  and  the  range  per  y-division  are 
computed  and  stored  in  positions  NDO  +  1  and  NDO  +  2  of  the  two  data  arrays. 
Throughout  this  portion  of  code,  NDO  is  the  number  of  time  series  points  (for 
each  file  )  to  be  plotted  on  the  present  page.  Headers  are  then  written, 
axes  formed,  and  the  present  page  plotted. 

The  user  is  now  asked  whether  he  wishes  a  Fourier  transform  performed, 
whether  he  simply  wishes  to  continue  with  the  time  series  plots,  or  whether 
he  wishes  to  abort  the  entire  process.  As  may  be  seen  from  the  listing, 
only  the  responses  'F1  and  'S'  cause  a  branch  to  code  that  has  not  already 
been  discussed.  In  the  case  of  the  'S'  response,  the  user  is  prompted  for 
two  variables ,  namely  PLTMIN  and  PLTSTP,  which  represent  the  minimum  value  of 
the  power  spectra  plot  (  in  decibels),  and  the  range  per  division.  After 
clearing  the  screen,  the  power  spectra  are  computed,  frequency  values  are  set 


« 


10 


up  in  XI,  headers  and  axes  are  plotted,  and  the  power  spectra  displayed.  We 
note  that  the  number  of  spectral  points  being  plotted  is  ND02,  which  may  not 
equal  one  half  of  NDO  since  the  time  domain  data  may  have  to  be  padded  to 
meet  the  fast  Fourier  transform  requirement  of  the  number  of  points  equalling 
a  power  of  two.  The  'F'  response  is  the  same  as  for  'S'  except  that  the 
user  is  not  prompted  for  PLTMIN  and  PLTSTP. 

The  program  then  pauses  to  permit  the  user  to  view  the  spectra,  and 
waits  for  a  response.  A  response  of  CR  causes  the  program  to  go  back  around 
the  main  plotting  loop,  whereas  a  response  of  'A'  causes  any  remaining 
plotting  to  be  ignored,  and  the  program  to  return  to  statement  number  884. 

Finally,  at  statement  number  998,  the  user  has  finished  with  the 
program,  the  initial  input  unit  is  closed,  any  temporary  files  created  by 
TEKPLOT  are  deleted,  and  the  code  exited. 

4.2.  SUBROUTINE  POWSPC 

This  routine  is  responsible  for  providing  a  fast  Fourier  transform, 
calculating  the  power  spectrum,  and  expressing  it  in  decibels. 

The  routine  first  finds  the  mean  from  all  of  the  data,  and  then  tapers 
the  data  with  a  sine  squared  function.  This  function  has  value  zero  at  the 

end  points  of  the  data,  and  has  increased  to  unity  over  a  range  of  10X  of  NDO 

within  each  end  of  the  data.  The  routine  then  finds  the  lowest  number  L 
which  is  a  power  of  two  and  greater  than  or  equal  to  the  number  of  data 
points.  The  data  of  length  NDO  is  then  padded  with  zeroes  until  it  is  of 
length  L. 

The  exact  call  to  the  fast  Fourier  transform  routine  is  system 

dependent  and  may  have  to  be  altered  for  different  installations.  In  our 

particular  system,  the  call  includes  the  area  from  which  data  are  to  be 
taken,  the  exponent  K  where  2**K  *  L,  an  array  S  into  which  the  internal  sine 
table  will  be  placed,  a  variable  indicating  that  the  sine  table  is  to  be 
calculated,  and  an  error  flag  IFERR. 
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On  return  from  the  transform  call,  error  conditions  are  checked  for, 
the  power  spectrum  computed,  and  the  result  expressed  in  dB.  The  routine 
then  exits. 

4.3.  SUBROUTINE  PLTOPEN 

This  routine  may  require  some  additional  lines  of  code  in  some 
installations.  In  its  present  form,  a  call  to  PLTNAME  initialises 

parameters  within  the  plot  package,  and  then  clears  the  screen.  The  call  to 
FINITT  forces  the  present  plot  buffer  to  be  output,  and  leaves  the  screen  in 
alphanumeric  mode  at  coordinate  (0,770).  Finally  the  call  to  limit  informs 
the  plot  package  of  the  maximum  logical  page  size. 

4.4.  SUBROUTINE  NWPLOT 

This  routine  has  entry  points  NWPLOT  and  NUPAGE.  Both  entries  clear 
the  screen  and  provide  a  suitable  time  lapse  before  returning  to  the  user. 
NWPLOT  first  pauses,  requesting  user  input,  before  clearing  the  screen.  A 
CR  response  causes  a  normal  return  to  the  calling  routine,  whereas  an  'A' 
response  causes  a  return  to  the  statement  indicated  in  the  calling  argument. 

4.5.  SUBROUTINES  HEADER  AND  FT AXES 

These  two  routines  contain  the  coding  necessary  to  put  most  of  the 
headers,  numbers  and  symbols,  on  both  the  time  series  and  power  spectra 
plots.  Except  for  AXISM,  all  subroutine  calls  are  standard  CALCOMP  calls. 
AXISM  is  a  common  modification  of  AXIS  which  suppresses  non-significant 
zeroes  from  the  axis  numerals  before  writing  them. 

The  entry  point  HEADR2  is  provided  to  enable  output  of  solely  the 
second  file  name. 


4.6.  SUBROUTINE  CLOSE 

This  subroutine  may  not  be  necessary  in  some  other  installations.  Its 

purpose  is  to  close  the  initial  input  unit,  and  delete  any  temporary  files 

that  may  have  been  created  by  TEKPLOT. 

4.7.  SUBROUTINES  BREAD  AND  BWRITE 

These  two  routines,  which  read  and  write  elements  of  a  data  file,  will 
almost  certainly  need  some  modification  to  satisfy  the  protocol  for  handling 
files  in  a  particular  installation. 

The  routine  BREAD  is  used  to  read  a  specified  number  of  elements, 

starting  at  a  particular  specified  point  in  a  file,  and  transfer  them  into  an 
array.  Within  a  particular  transfer,  the  elements  are  sequential,  but  a 
subsequent  call  to  BREAD  might  ask  for  the  transfer  of  a  block  from  a 

completely  different  part  of  the  file.  The  routine  handles  reads  from  both 
tape  and  disk.  In  the  present  implementation,  the  routine  is  structured  to 
read  from  permanent  disk  files  (ITAPE(I)  *  9998),  or  temporary  disk  files 
(ITAPE(I)  =  9999). 

The  routine  BWRITE  essentially  performs  the  inverse  function  of  BREAD, 
transferring  from  an  array  into  an  output  file.  The  coding  however,  need 
only  support  write  operations  to  the  temporary  disk  (9999),  since  the 
subroutine  is  never  called  upon  to  write  to  permanent  disk. 

Any  system  dependent  calls,  needed  to  support  tape  labels,  file 
formats,  file  searching  etc.,  should  be  placed  in  these  two  routines. 


APPENDIX  A 


1.  Flow  Chart  for  TEKPLOT  Main  Routine 
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Read  data  for  one  page  full 
of  the  time  series. 
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Flow  Chart  for  Power  Spectrum  Routine 
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TKKPLO? 

P.M.  HOLT  HA.” 

D.R.E.P. 

VICTORIA,  R.C. 
CANADA 

NOVEMBER  la0! 


TKKPLOT  ENABLES  INTERACTIVE  VIEWING  OF  UP  TO  TWO 
FILES  CONTAINING  TIME  SERIES  DATA.  FILES  CAM  RFSIDE 
EITHER  ON  MAGNETIC  TAPE  OR  ON  DISK. 

THE  USER  CAN  VIEW  SELECTED  PORTIONS  OF  HIS  DATA,  AND 
AND  CAN  CHOOSE  FURTHER  SEGMENTS  WHICH  ARE 
EITHER  FURTHER  AHEAD,  OR  BEHIND,  IN  THE  FILE.  SEVERAL 
SCALING  MODES  ARE  AVAILABLE  FOR  PLOTTING  THE  DATA 

FOURIER  TRANSFORMS  OF  THE  SEGMENTS  CAN  BE  CALCULATED , 
AND  THE  RESULTANT  TOWER  SPECTRA  PLOTTED. 


REAL  DATA(  U098 )  ,X1  ( U098 )  ,DAT2  ( 1*098 ) 

COMMON  ND02 , I DUM , DATA , X 1 , D AT2 
LOGICAL  DOTIME , TWO PI.T , DISK 
DIMENSION  NTAP2(5),TEMPDAT(U) 

INTEGER  NTAPE(5) ,IFILE(S) ,IFILE2(5) 

COMMON  /HDR/  NTAPE,IFILE,TFILE2,NST,IH,IM, IS, NSEC, DOTIME 
COMMON  /FTCOMM/  PLTMI N , PLTSTP , AXS , I SP ,TVOPLT 
DATA  NTAP2 ( 1 ) /9999 / 

INTEGER  A/ ' A' / ,FEFF / ' F’ / ,ISPAC / *  '/ 

DATA  RMIN , RMAX/Z7FFFEFFF, 7.80000001/ 


SET  UP  MAXIMUM  NUMBER  OF  POINT’S  PER  PAGE 
NMAX=U09h 


OPEN  THE  PLOTTER  AND  CLEAR  THE  SCREEN 
CALL  PLTOPEN 


PROMPT  THE  USER  FOR  INPUT  PARAMETERS  AND  READ  THEM 
PARAMETERS  ARE: 

FIRST  KILE  NAME 
TAPE  OR  DISC  NUMBER 
SECOND  FILE  NAME 
START  POINT  IN  THE  FILE(S) 
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52. 

r 

TOTAL  NUMRF.R  OF  SECONDS  OF  DATA  TO  BE  PLOTTED 

5  3. 

■’ 

NUMBER  OF  SECONDS  TO  BF.  PLOTTED  ON  EACH  PACE 

5>*. 

r 

SAMPLE  RATE  OF  DATA  ( INTEGER) 

SB. 

0 

STARTIN'".  TIME  OF  THE  RFC. INN I NO  OF  FILE(s) 

y>. 

r 

SCALING  MODE 

57. 

c 

SCALE  PARAMETERS  IF  SCAJ.E  MODE  EQUALS  2 

5B. 

c 

SQ. 

WRITE(6 ,101 ) 

60. 

101 

FORMAT (  '  ENTER  FILE  NAME  IN  2A4 ' ) 

61. 

READ( 5 ,83  3  ( I?ILE( I ) ,1=1 ,2) 

62. 

83 

FORMAT ( 2 AU) 

63. 

WRITE(6,292) 

64. 

202 

FORMAT(/'  ENTER  TAPE  ,  14,  3A4 ,  CR  DEFAULTS  TO  9998’) 

65. 

READ(5,293)  ( NTAPE(  J ) , J=1 ,4  ) 

66. 

293 

FORMAT ( I 4, 3 A4 ) 

67. 

WRITE(6,202) 

68. 

202 

FORMAT(/'  FILF  2.  CR  =  NONE’) 

69. 

READ(6,83)  ( IFILE2(l) ,1=1,2) 

70. 

884 

WRITE(6,Bll) 

71. 

811 

FORMAT(/'  ENTER  NSTART,  TIME  LENGTH,  AND  SECS  PER  PAGE’, 

72. 

1  r  DEFAULTS  ARE  1,  20,  AND  10',’  FORMAT  110') 

73. 

READ  ( 5 ,80 ,END=903 )  NSTART ,NDISP , NSEC 

74. 

80 

FORMAT (110) 

75. 

write(  6,291+) 

76. 

29  U 

FORMAT ( / '  ENTER  SAMPLE  RATE  IN  15,  CR  IS  40  SAMP  PER’, 

77. 

1  ’  SEC') 

78. 

READ(5 ,295 )  ISP 

79. 

296 

FORMAT (15) 

80. 

WRITE(6,771) 

81. 

771 

FORMAT ( / '  ENTER  START  TIME  OF  TAPE.  CR  TO  IGNORE', 

32. 

1  '  THIS  OPTION.  12,213*) 

83. 

READ(5 ,772)  I HST , IMS  T , ISST 

84. 

772 

FORMAT( 12,213) 

89. 

TH 

WRITE( 6,77U) 

86. 

774 

FORMAT ( / '  ENTER  SCALING.  CR  FOR  PER  PAGE,' 

87. 

1  '  1  FOR  WHOLE  DISPLAY,  2  TO  BF  INPUT’) 

88. 

READ( 5 ,775 )  ISCALE 

80. 

775 

FORMAT (II) 

90. 

n 

LOOP  BACK  IF  ILLEGAL  SCALE  MODE 

91. 

IF( ( ISCALE. EQ.O ) .OR. ( ISCALE. EQ. I ) .OR. ( ISCALE. EQ.?)  V.OTO  877 

92. 

GOTO  777 

93. 

377 

CONTINUE 

94. 

IF  ( ISCALE. EQ.l)  WRITE( 6 ,778 ) 

95. 

778 

FORMAT ( /'  SCALE  THROUGH  TOE  WHOLE  DATA  ’) 

06. 

IF  (ISCALE. EQ.O)  WRITE(6,77o) 

97. 

770 

FORMAT ( / '  SCALE  PEP  PAGE  ’ ) 

98. 

c 

99. 

IE  ( ISCALE. NE. 2)  GOTO  878 

100. 

r* 

ENTER  SCALING  PARAMETERS  IF  ISCALE=2 

101. 

write(6 ,879) 

102. 

870 

FORMAT(/'  ENTER  Y  SCALE,  !JNITS  PER  INCH,  F10.5’) 

103. 

READ(5 ,890)  YS 

io4. 

800 

FORMAT( F10.5 ) 

105. 

V/RITE(6 ,891 )  YS 
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106. 

891 

107. 

108. 

878 

109. 

110. 

r 

111. 

C 

112. 

1’ 

113. 

C 

llU. 

115. 

116. 

117. 

118. 

119. 

120. 

121. 

L22. 

c 

123. 

l?.k. 

125. 

126. 

127. 

P 

128. 

129* 

c 

130. 

c 

131. 

c 

132. 

r 

133. 

13U. 

135. 

C 

136. 

137. 

138. 

139. 

101 

1**0. 

lUl. 

1**2. 

c 

1U3. 

100 

lUL. 

c 

1*»5. 

1U6. 

11*7. 

c 

1U8. 

1U9. 

150. 

c 

151. 

c 

152. 

c 

153. 

151*. 

81» 

155. 

156. 

203 

157. 

158. 

773 

159. 

FORMAT ( / '  Y  SCALE  WILL  PE  •  ,01:'. 5,'  UNITS  PER  INCH',/ 

1  Y  MIN  WILL  BE  CHOSEN  AUTOMATICALLY  FOR  EACH  GRAPH' ) 
CONTINUE 


SET  UP  DEFAULT,  CONTROL  AND  SUNDRY  OTHER  PARAMETERS 
IF  (NTAPE(l).EQ.O)  NTAPE; 1 )=999B 
DISK*. FALSE. 

IF  ( (NTAPE( 1 ). EQ. 9999). OR. (NTAPE(l).EQ. 999*0)  DISK*. TRUE. 
TWOPLT* . TRUE . 

IF  ( IFILE2( 1 ) . EQ. ISPAC )  TWOPLT*. FALSE. 

IF  (ISP.EQ.O)  ISP=LO 
DOTIME*. TRUE. 

IF  ( I HST+ IMST+ I SST . EQ . 0 )  DOTIME*. FALSE. 

NPNTS  IS  TOTAL  NUMBER  OF  POINTS.  'TO  BE  DISPLAYED 
NPNTS=NDISP*ISP 
IF  ( NSTART.EQ.O )  N START* 1 
IF  ( NPNTS. EQ.O)  NPNTS=800 
IF  (NSEC. EQ.O)  NSEC=10 

ISTART  IS  START  TIME  OF  TAPE  EXPRESSED  IN  SECONDS 
ISTART=60*(IHST#60+IMST)+ISST 
PAGE  CONTROL  PARAMETERS 
NOIN  IS  NUMBER  OF  INCHES  PER  PAGE 
N POINT  IS  NUMBER  OF  POINTS  TO  BE  PLOTTED  EACH  PAGE 
PERIN  IS  THE  POINT  DENSITY  PER  INCH 
NOIN=10 

NP0INT=NSEC*I3P 
CHECK  N POINT  NOT  TOO  LARGE 
IF  ( NPOINT.LE.NMAX)  GOTO  100 
CALL  NUPAGE 
WRITE(6,10l)  NPOINT 

FORMAT(/  '  NPOINT  OF  ',120,'  "00  LARGE.  RUN  ABORTED') 

CALL  NWPLOT( B8LS) 

GOTO  88 L 

PERIN=NPOINT*0 . 1 

SET  UP  DEFAULT  POWER  SPECTRA  PLOT  PARAMEYERS 

pltmin=-6o. 

PLTSTP=10. 

SET  UP  AXS  EQUAL  TO  NUMBER  OP  PLOTS  TO  BE  DONE 
AXS=1 

IF  (TWOPLT)  AXS=2 

OUTPUT  PARAMETERS  FOR  USER  REFLECTION 

WRITE(6,84)  (NTAPE(I),I*1,1),IFII.K( l),IFTLE(2) 

FORMAT(  '  TAPE  '  ,!’>  FILE  '  ,2AM 

IP  (TWOPLT)  WRITE(6,20  5 )  (  IFILE?(  I )  ,1=1  ,.?) 

FORMAT( 21X ,2AU ) 

IF  (DOTIME)  WRITE  (6,773)  IHGT.IMST.IGST 
FORMAT (  *  TAPE  START  ’ ,3X,3I3) 

WRITE(6 ,88 )  NSTART , NPNTS ,NSEC , PERI N 


mntmam 


m 
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160. 

32 

FiM\MAT(  '  RTAF’T  AT  '  ,5X,T8,«  NUMBER  OF  !>0INT9>  '  ,T8, 

Ihl  . 

1  /'  SECONDS  PER  PAGE  ',TS,'  DATA  POINTS  PER  TNCH 

169. 

9  F6 . 2 ) 

16 

0 

lnh. 

0 

ISS. 

r 

1  tS6 . 

U 

It  >7. 

r 

168. 

C. 

IF  FILE  IS  NOT  ALREADY  ON  DISK,  OOPY  FI LEU  FROM 

l6o. 

c 

TAPE  TO  TEMPORARY  DISK,  TO  AVOID 

170. 

c 

MULTIPLE  TAPE  REWINDS 

171. 

c 

172. 

,7=0 

173. 

IF  ( .NOT.TWOPLT)  OOTO  39? 

17U. 

IF  (DISK)  OOTO  393 

175. 

IF  (NPNTS.LT. L0q6)  OOTO  361 

176. 

p 

TRANSFER  IN  RLOCKS  OF  L096 

177. 

DO  353  I=1,NPNTS,L0Q6 

178. 

IF  (I+LO96-I.OT.NPNTS)  GOTO  361 

179- 

J=J+1 

l8n. 

CALL  BREAD ( NTAPE , I FILE2 , L096 , DAT2 , I- 1+NSTART ) 

lfll. 

CALL  BWRITE(  NTAP2 , 1 FTLE2 ,L096  ,DAT2 , 1-l+NSTART ) 

182. 

353 

CONTINUE 

163. 

361 

CONTINUE 

m. 

NUM=NPNTS-J*L006 

185. 

IF  (NUM.EQ.O)  OOTO  398 

136. 

0 

TRANSFER  ANY  REMAINDER 

187. 

CALL  RREAD( NTAPE , IFILE2 ,NUM,DAT? , J*L096+NSTART ) 

133. 

CALL  RWRITE ( NTAP2 , IFILE2 , NUM , DAT2 , J*L096+NSTART ) 

139. 

398 

CONTINUE 

190. 

n 

191. 

r 

192. 

c 

193. 

c 

19U. 

c 

HANDLE  SCALE  MODE  =  1  SITUATION 

195- 

0 

LOOK  FOR  MAXIMA  AND  MINIMUM  VALUES  OF  THE  SPECIFIED 

196. 

c 

TIME  SERIES  FOR  FILE  ONE 

107. 

IF  ( ISCALE.NE. 1 )  GOTO  8oi 

198. 

J=0 

109. 

IF  (NPNTS.LT.L006)  GOTO  261 

900. 

c 

SCAN  THROUGH  DATA  IN  RLOCKS  OF  LO96 

201. 

DO  253  I=1,NPNTS,L096 

902. 

IF  (  I+L096-1  .GT.NPNTS  )  GOTO  26l 

203. 

J=J+1 

90  U. 

CALL  BREAD( NTAPE, IFILE,L006, DATA, I-l+NSTART ) 

905. 

0 

COPY  FILE  ONE  TO  DISK  IF  IT  WAS  NOT  ALREADY  THERE 

906. 

IF  (.NOT. DIRK)  CALL  RWRITE( NTAP2 ,IFILE,L096 ,DATA , I-l+NSTART) 

207 . 

DO  257  K=1,L096 

208. 

IF  (DATA(K).GT.RMAX)  RMAX=DATA(:<) 

209. 

IF  (DATA(K).LT.RMIN)  RMIfJ=DATA( K ) 

910. 

257 

DATA(K)=1.E35 

211. 

253 

CONTINUE 

212. 

261 

CONTINUE 

213. 

NUM=NPNTS-J«L0o6 

214. 

c 

SCAN  THROUGH  ANY  REMAINING  POINTS 

215. 

IF  (NUM.EQ.O)  GOTO  108 

216. 

CALL  BREAD ( NTAPE , I FI LE ,NUM .DATA , J*4096+NSTART ) 

217. 

IF(. NOT. DISK)  CALL  BWRITE(NTAP2,IFILE,NUM,DATA,J*4096+NSTART) 

213. 

do  259  k=i,num 

219. 

IF  ( data(k) .GT.RMAX )  RMAX=DATA(K) 

220. 

IF  (DATA(K).LT.RMIN)  rmin=data(k) 

221. 

259 

DATA(k)=1.E35 

222. 

198 

CONTINUE 

223. 

c 

CALCULATE  THE  MINIMUM,  AND  THE  RANGE  PER  DIVISION 

224. 

c 

NEEDED  FOR  THE  PLOT 

225. 

TEMPDAT(1)=RMAX 

226. 

TEMPDAT(2)=RMIN 

227. 

CALL  SCALE  (TEMPDAT,3. ,2,l) 

228. 

YMIN=TEMPDAT(  3 ) 

229. 

YSEP=TEMPDAT(4) 

230. 

WRITE(6,1974)  RMIN,RMAX 

231. 

1974 

FORMAT(  '  MIN  AND  MAX  SIGNAL  ’ ,2E15.5) 

232. 

WRITE(6,150)  ymin,ysep 

233. 

150 

FORMAT(  '  MINIMUM  Y  AND  SEPARATION  ' ,2E15.5) 

234. 

301 

CONTINUE 

235. 

0 

236. 

0 

237. 

c 

233. 

c 

239- 

c 

240. 

c 

BEGINNING  OF  THE  DATA  PLOTTING  LOOP 

241. 

c 

TWO  ENTRY  POINTS:  10  AND  110,  DEPENDING  ON  WHETHER 

242. 

c 

CLEAR  PAGE  NEEDS  A  PAUSE  OR  NOT 

243. 

START=NSTART 

244. 

NDO=N POINT 

245. 

T0TD0N=0 

246. 

NST=NSTART-N POINT 

247. 

J=0 

248. 

10 

CALL  NWPL0T(884S) 

249. 

GOTO  111 

250. 

110 

CALL  NUPAGE 

251. 

c 

INCREASE  NUMBER  OF  POINTS  DONE  BY  A  PAGE 

252. 

111 

TOTDON=TOTDON+N PO I NT 

253. 

c 

IF  LAST  PAGE  PLOTTED  WAS  NOT  Ft  ILL,  THEN  MUST  HAVE  FINISHED 

254. 

c 

THE  TIME  SERIES.  GO  BACK  AND  ASK  FOR  ANOTHER 

255. 

c 

TIME  SERIES  SEGMENT 

256. 

IF  (NDO.LT.NPOINT)  GOTO  884 

257. 

c 

INCREASE  STARTING  POINT  BY  ONE  PAGE 

253. 

NST=NST+N POINT 

259. 

c 

IF  A  COMPLETE  PAGE  WOULD  TAKE  US  PAST  THE  END  OF 

260. 

c 

THE  SPECIFIED  TIME  SERIES,  ADJUST  NUMBER  OF  POINTS 

261. 

IF  (TOTDON.GT.NPNTS)  NDO=NPOINT-(TOTDON-NPNTS) 

262. 

c 

IF  NDO  IS  ZERO,  THEN  THE  WHOLE  TIME  SERIES  SEGMENT 

263. 

n 

HAS  BEEN  PLOTTED  EXACTLY,  WITHOUT  THE  NEED  FOR  A 

264. 

c 

SHORT  PAGE.  GO  BACK  AND  ASK  TOR  A  NEW  TIME  SERIES 

265. 

0 

SEGMENT 

266. 

IF  (NDO.EQ.O)  GOTO  884 

267. 

0 

INCREASE  POINTER  TO  END  OF  DATA  IN  CURRENT  PLOT 

A9 


;>6R. 

P6'».  c 

To.  c 

071  . 

■  >•  S 
I  ■  • 

OT  l. 

'  r’i . 

?7r,.  (> 

vro.  c 

OTT. 

07,9. 

070. 

:’Ro.  Rl'i 

281.  o 

282.  n 

28  3. 

0BI4 . 

285. 

-86. 

?q7.  20 

2?8. 

289. 

299 .  9 

291.  0 

292.  0 

293.  ('  SET  UP  THE  SCALING  APPROPRIATE  TO  THE  THREE  MODES 

29  L .  C 

29r).  IE  ( ISC  ALE.  NE.  1 )  GOTO  804 

296.  0  MODE  1 

297.  0  ALL  PAGES  OF  ROTH  FILES  ARE  SCALED  IDENTICALLY. 

998.  C  SCALING  IS  APPROPRIATE  TO  THE  MAXIMUM  AND  MINIMUM  VALUES 

299.  C  OF  THE  WHOLE  OF  THE  SPECIFIED  TIME  SERIES  OF  FILE  ONE 

300.  DATA(NDO+l)=YMIN 

301.  DATA ( NDO+2 )=YSEP 

302.  DAT2 ( NDO+ 1 ) =YMI N 

303.  DAT2 ( NDO+2 ) =YSSP 

30U.  GOTO  806 

305.  C 

306.  C 

307.  80L  IF  ( ISCALE.NE.O )  GOTO  805 

308.  C  MODE  0 

309.  C  EVERY  PAGE  IS  SCALED  SEPARATELY  AND  AUTOMATICALLY. 

310.  C  BOTH  FILES  ARE  SCALED  INDEPENDENTLY 

311.  CALL  SCALE( DATA ,8. ,ND0,1 ) 

312.  IF(TWOPLT)  CALL  SC4LE( DAT2,R. ,ND0,1 ) 

313.  GOTO  806 

31*4.  C 

315.  C 

316.  C  MODE  2 

317.  C  RANGE  PER  DIVISION  WAS  READ  IN,  AND  IS  USED  THROUGHOUT 

318.  C  MINIMA  ARE  FOUND  SEPERATELY  FOR  EACH  PAGE  AND  FOR  EACH  FILE 

319.  805  CALL  SCALE(DATA,8. ,ND0,1) 

320.  IF  (TWOPLT)  CALL  SCALE( DAT2 ,8 . ,NDO ,1 ) 

321.  DATA( NDO+2 )=YS 


NFLN=NST+NP0-1 

HEAD  THE  DATA  FOR  FILE  ONE 
TP  ( ( . NOT. DISK'. AND. (I3CALE. EQ. 1) ) 

I  CALL  BREAD(NTAP2,IFILE,NDO,DATA,NST,0) 

IF  ( (DISK). OR. ( ( .NOT. DISK). AND. ( ISCALE. NE. 1 ) ) ) 

1  CALL  RREAP(Nr’APE,TEILR,NDO,D£TA,NST) 

READ  DATA  FOR  SECOND  FILE  IF  NEEDED 
IF  (.NOT. TWOPLT)  GOTO  RlL 

TF( .NOT. DISK)  ''ALL  BREAD ( NTAP2 , TFILE2 ,NDO ,DAT2 , NST ,0 ) 
IF( DISK'  CALL  :SREAD(NTAPE,IFILE2,ND0,DAT2,NST) 
CONTINUE 

SET  UP  X-AXIS  POINTS  AND  THE  SCALE  INFORMATION 
1=0 

DO  20  IK=NST,NFIN 
1=1+1 
Xl( I)=IK 
CONTINUE 

XI  ( NDO+1 )  =FLOA'n(  NST ) 

XI ( NDO+2 )=PERIN 
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322. 

DAT2 ( NDO+2 ) =Y3 

323. 

C 

32k. 

P 

325. 

P 

HALVE  THE  SCALE  IF  TWO  PLOmS  TO  BE  DONE 

326. 

Bo  6 

IF  (TWOPLT)  DATA (NDO+2) =DATA ( N  DO+2 ) *2 . 

327. 

IF  (TWOPLT)  DAT2 ( NDO+2 )=DAT2( NDO+2) *2. 

328. 

329. 

0 

CALCULATE  TIME  AT  START  OF  CURRENT  PAGE 

330. 

NSECDON=  (NST-1 ) /ISP  +  I3TART 

331. 

ih=nsecdon/36oo 

332. 

tt,eft=nsecdon-ih»36oo 

333. 

im=ileft/6o 

33**. 

is=ileft-im*6o 

335. 

IF  (lH.GE.2l*)  IH=IH-2l* 

336. 

r 

, 

337. 

0 

WRITE  HEADER  FOR  THE  TIME  SERIES  PET'1 

338. 

CALL  HEADER (TWOPLT) 

339. 

CALL  SYMB0L(7.50 ,7.0, .IS , ’SCALE  MODE’ , .0,11) 

3 1*0. 

CALL  NUMBER  (9. 5 ,7.0  ,  .15  ,  FLOAT  (  [SCALP;  ),.0,0) 

3**1. 

c 

3l*2. 

0 

3l*3. 

c 

DO  THE  TIME  SERIES  PLOTS 

3UU. 

IF  (TWOPLT)  CALL  PL0T( .0 ,U.  ,-3 ) 

31*5. 

c 

DRAW  THE  AXIS  FOR  THE  FIRST  FILE  DATA  PLOT 

31*6. 

c 

EIGHT  INCHES  IF  ONE  PLOT,  THREE  INCHES  IF  TWO  PLOTS 

3U7. 

AXL=8. 

31*8. 

IF  (TWOPLT)  AXL=3- 

349. 

CALL  AXISM( .0 , ,0 , ’  ’ ,1 ,AXL,90 .  ,DATA(ND0+1 ) ,DATA( NDO+2) 

350. 

1  ,.15, .15) 

351. 

c. 

PLOT  THE  DATA 

352. 

CALL  LINE  ( XI, DATA, NDO, 1,0,0) 

353. 

CALL  FINITT(0 ,0 ) 

351*. 

c 

355. 

cc 

356. 

c 

DO  THE  SECOND  FILE  IF  REQUIRED 

357. 

IF  (.NOT. TWOPLT)  GOTO  11 

358. 

CALL  PLOT(  .0,-1*.  ,-3 ) 

359. 

c 

DRAW  THE  AXIS  FOR  THE  SECOND  DATA  PLOT 

360. 

CALL  AXISM( .0 , .0 , ’  * ,1 ,AXL ,90 . ,DAT2( NDO+1 ) ,DAT2( NDO+2) 

36l. 

1,  .15  15) 

362. 

c 

PLOT  THE  DATA 

363. 

CALL  LINE( XI  ,DAT2, NDO, 1,0,0) 

36U. 

11 

CONTINUE 

365. 

CALL  FINITT(0 ,0 ) 

366. 

c 

367. 

c 

368. 

c 

THAT’S  THE  TIME  SERIES  PLOTTING  FINISHED,  WHAT  NEXT. 

369. 

c 

370. 

c 

FOURIER  TRANSFORM  REQUIRED? 

371. 

WRITF.(6 ,12) 

372. 

L2 

FORMAT ( ’  ENTER  CR  FOR  CONTINUE,  F  EOR  F.T.,  S  YOU  SCALE ’ , 

373. 

1  * ,  A  FOR  ABORT’ ) 

37U. 

CALL  FINITT(860,0) 

375. 

READ( 5,13  ,END=998 )  IRESP 

All 


<  ’n. 

377. 

1 3 

FORMAT (  A1 ) 

7 

A  =  10*100 K  CPU, -UN  I  NS  OF  PLOTS,  OR  =  IGNORE  F.F.T. 

V’O. 

f  ok  s  1: -pi. ,v-  a  "ransform 

;Ho. 

1  v  ( TRESr.TO. 0  TAIL  T'PAGF 

3*1. 

IF  (  I'?FSP.  EO.  a)  10TO  R°.L 

<8  1 . 

if  ( rHFFp.Fii.rr'PAO)  goto  no 

;•??. 

?'V.. 

?3s. 

(P  ' [PROP, EC. f^FF'  9070  1012 

i' 

0 

l’,;V’R  VI.  1  >7  SUPPLY  POWER  Appf'TRA  PLOT  SCALING 

ClA. 

FALL  T'PASF 

337. 

WRITF(6 ,101? I 

3R9 . 

101? 

FORMA'" f '  ENTER  PI  '7  MINIMUM  IN  DB.  F10.5') 

l^Q. 

READ{  5 ,101>4 )  PLTMIN 

<00 . 

10 1. 1 

FORMAT ( F10 .5 I 

301. 

WRITER  >1915 ) 

3<??. 

1015 

FORMAT ( ’  ENTER  PUT  SCALE  PER  INCH  IN  OB.  F10.5') 

393. 

READ(  5  ,101>< '  PLTSTP 

3Q1* . 

395. 

n 

CALCULATE  7!!F  POWER  SPF/'TR' JM  AFTER  CLEARING  THE  SCREEN 

yyfu 

1012 

CALI.  NO  PAGE 

307. 

CALL  POWSPr  (  ?  A  “"A  ,NPO ) 

7)ft. 

r* 

PUT  Ill  THE  SCAT  INC 

399. 

DATA( ND02+1 )=PLTUIN 

■'<90. 

oAmAf  N  1)0  2 + 2 }  =  PLTS'  I AX  3 

>401 . 

r 

>402. 

0 

U03. 

c 

PUT  HEADER  AND  AXES  ON  F.T.  PLOTS 

'»0H. 

r 

FORCE  THE  SECOND  FILE  NAME  TO  BE  OMITTED 

HOS . 

CALL  HEADER(. FALSE.) 

•»0'5  • 

CALL  FTAXES 

U07 . 

r 

PTO  ON  ANY  SECOND  FILENAME  AFTER  PLOT  HAS  BEEN 

'.09. 

c 

FACTORED  BY  FTAXES  CALL 

U09. 

IF  ( TWOPLT )  CALL  HFADR2 

HlO. 

c 

•'.11 . 

c 

Hi?. 

c 

SET  UP  THE  X-AXTS  AND  SCALING 

413. 

DO  703  IX=1,ND02 

J»1  U . 

703 

X1(IK)  =  .5*FL9AT(ISP)#FL0AT(IK-1)/'FL0AT(ND02) 

U15. 

Xl(?ID02+l)=0. 

M6. 

Xl(ND02+2  )=.f'5*ISP 

Ul7. 

C 

'il3. 

C 

PLOT  THE  FIRST  POWER  SPECTRUM 

'4 19. 

IF  (TWOPL"1}  CALL  PLOTf  .9,U.5,-3) 

U20. 

CALL  LINE( XI .DATA ,ND02 ,1 ,0 ,0 ) 

U21. 

c 

>i22. 

7 

CALCULATE  AND  PLOT  THE  SECOND  FILE  SPECTRUM  IF  NEEDED 

U23. 

IF  (.NOT. TV/O PLT)  GOTO  hi 

liO'i. 

CALL  ROWS PC ( D AT2 , NDO ) 

U25. 

DAT2 ( ND02+1 ) =PLTMIN 

U26. 

9AT2( ND02+2 }= PLTSTP* A XS 

I'27  * 

CALL  FINITT(30,30’' 

U2B . 

CALL  PLOT(.G,-Ii.5,-3) 

>4  29. 

CALL  LINE( XI, DAT2,ND02, 1,0, 0) 

A12 


1*30. 

r 

1*31. 

r 

1*32. 

c 

POWER  SPECTRA  COMPLETED.  00  BACK  AROUND  THE  PACE  LOO? 

U33. 

hi 

CALL  FI MITT (30, 30) 

U31*  • 

GOTO  10 

1*35. 

C, 

1*36. 

1*37. 

a 

USER  HAS  FINISHED  WITH  THE  PROGRAM 

1*38. 

ij 

CLOSE  INPUT  UNIT  AMD  DELETE  AMY  TEMPORARY  FILES 

1*39. 

c 

CREATED  BY  TEKPLOT 

1*1*0. 

098 

CALL  CLOSE( NTAPE,NTAP2 ,DISK) 

1*1*1. 

1^ 

1*1*2. 

STOP 

1*1*3. 

EMD 

1. 

SUBROUTINE  P0WSPC( A ,NDO) 

2. 

DIMENSION  A(1),S(2056) 

3. 

COMMON  ND02 

1*. 

COMMON  /FTCOMM/  PLTMIN,PLTSTP,AXS,ISP,TWOPLT 

5. 

C 

6. 

C 

7. 

c 

8. 

c 

SUBTRACT  OUT  DC  AFTER  FINDING  THE  MEAN 

9. 

SUM=0 . 

10. 

DO  3  IK=1 ,NDO 

11. 

SUM=SUM+A( IK ) 

12. 

3 

CONTINUE 

13. 

SUM=SUM/NDO 

lU. 

DO  U  IK=1 ,NDO 

15. 

4 

A(IK)=A( IK ) -SUM 

16. 

C 

17. 

c 

WINDOW  DATA  OVER  10*,  EACH  END  WITH  SIN*»2  WINDOW 

18. 

NWIND=ND0/10  +  1 

19. 

c 

SKIP  IF  NO  POINTS  TO  WINDOW 

20. 

IF  (NWIND.EQ.O)  GOTO  5555 

21. 

DO  255  IS=1,NWIND 

22. 

rsin=sin( . 5*3. i4 15927* ( IS-l) /(NWIND-l) ) 

23. 

RSIN=RSIN*RSIN 

21*. 

a( is )=a( is)*rsin 

25. 

255 

A(ND0+1-IS)=A(ND0+1-IS)*RSIN 

2  6. 

c 

27. 

C 

FIND  POWER  OF  2  GREATER  OR  EQUAL  TO  NDO 

28. 

5555 

L=1 

29. 

DO  1  1=1,20000 

30. 

K=I 

31. 

L=2*L 

32. 

IF  (L.GE.NDO)  GOTO  2 

33. 

1 

CONTINUE 

3U. 

C 

L  IS  AN  EXACT  POWER  OF  2,  AND  IS  EITHER  EQUAL  TO  NDO 

35. 

C 

OR  THE  NEXT  POWER  OF  2  ABOVE  NDO 

36. 

C 

o  o  o  -^i  noon 


<'  PAD  WITH  ZEROES  I?  L  WAS  NOT  EQUAL  TO  NDO 

2  NUMZ=L-NDO 

IP  (NUMZ.EQ.O)  goto  5S 

ndoi=ndo+i 

no  5  IK=ND01,L 
5  A( IK )=0. 

c 

CALCULATE  THE  FFT  OF  ARRAY  A,  2**K  ELEMENTS  EXACTLY 
C  S  IS  STORAGE  FOR  THE  SINE  TABLE 

C  -1  SAYS  CACULATE  mH£  SINE  TABLE  EACH  TIME 

C  IFERR  IS  AN  ERROR  FLAT.  0  =  OK 

c  -1  =  SINE  TABLE  ALREADY  DONE 

55  CALL  CORRFFTf  A ,K ,S ,-l , IFERR ) 

n 

C  CHECK  FOR  ANY  ERROR  CONDITION 

I p  ( -NOT. ((IFERR. EQ.O). OR. (IFERR. EQ.-l)))  WRITE(6,6)  IFERR 
FORMAT ( / 1  ERROR  12 , '  IN  F.T.  »*»*»*»#»***##') 

SET  UP  NUMBER  OF  FOURIER  COEFFICIENTS  TO  PLOT 
NDO 2=  L/2  +  1 

CALCULATE  THE  POWER,  AND  BOTTOM  JUSTIFY  IT  IN  ARRAY  A 
A(l)=A(l)**2 
TEMP=A(2)**2 
DO  7  IK=3,L,2 

A( (lK+l)/2)=  (A(lK)**2+A(lK+l)**2) 

continue 

A(ND02)=TEMP 


FORM  POWER  IN  DB 
DO  8  IK=1  ,ND02 

a( ik)=io.*alogio(a(ik) ) 

c 

B  IK  (A(IK).LT.PLTMIN)  A(IK)=PLTMIN 
RETURN 
END 


SI JB ROUTINE  PLTOPEN 
C 

C  THIS  SUBROUTINE  INITIALIZES  PLOTTING  VARIABLES 
C  AND  CLEARS  THE  SCREEN 

r* 

CALL  PLTNAME( 'MYNAME' ,6) 

CALL  FINITT(0,7?0) 

CALL  LIMIT(0. ,10. ,0. ,lU.) 

RETURN 

END 
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SUBROUTINE  NWPLOT**) 

INTEGER  A/' A'/ 

INTEGER  ISPAC/*  '/ 

CALL  FINITT(0 ,0 ) 

WRITE(6,10) 

id  formatC  enter  cr  to  continue,  a  to  abort') 

CALL  FINITT* 480,0) 

READ* 5, 12)  ICR 
12  FORMAT(Al) 

GOTO  20 

ENTRY  NUPAGE**) 

ICR=ISPAC 
20  CONTINUE 

CALL  FINITT{ 0,48) 

CALL  NEW PAG 
CALL  TSEND 
CALL  PLTNAME(O.O) 

CALL  LIMIT*0.  ,10.  ,.0,l4.) 

IF  (ICR.EQ.A)  RETURN  1 
RETURN 
END 


I 


[ 


SUBROUTINE  HEADER (TWOPLT) 

DIMENSION  NTAPE{ 5 ) , IFILE( 5 ) .IFILE2* 5 ) 

LOGICAL  TWOPLT, DOTIME 

COMMON  /HDR /  NTAPE,IFILE,IFILE2 ,NST,IH ,IM, IS .NSEC , DOTIME 
CALL  SYMBOL* 4. 5, 7. 5,. 15,  ’FILE  ’,.0,5) 

CALL  SYMBOL(6. ,7.5,.15,IFILE(l),.0,3) 

CALL  SYMBOL  (7.5 ,7. 5 , .15 , 'TAPE  ’,.0,5) 

CALL  NUMBER  (8.5 ,7.5 , . 15 ,FLOAT(NTAPE( 1 ) ) , .0 ,0 ) 

CALL  SYMBOL^. 5, 7. 25, .15,  ' NSTART  ’,.0,8) 

CALL  NUMBERS. , 7. 25, .15, FLOAT* NST),. 0,-1) 

IF  (.NOT. DOTIME)  GOTO  800 
CALL  SYMBOL  (7.5 ,7.25 ,.15 , ’TIME  ’,.0,5) 

CALL  NUMBER  (8.5,7.25,.15,FLOAT(IH) ,.0,-l) 

CALL  NUMBER  (9* ,7*25 , .15 ,FL0AT(IM) , .0 ,-l) 

CALL  NUMBER  (9*50 ,7.25 ,. 15 .FLOAT* IS) ,.0,-1) 

800  CONTINUE 

CALL  SYMBOL*^. 5, 7.0,. 15, 'NSEC  ’,.0,5) 

CALL  NUMBER**).  ,7.  ,.15  .FLOAT* NSEC)  , .0,-1) 

C 

C  IF  ONLY  ONE  PLOT  TO  DO  THEN  RETURN 
C  ELSE  WRITE  THE  SECOND  FILE  NAME 
C 

IF  (.NOT. TWOPLT)  RETURN 
C 

C  ENTRY  POINT  FOR  JUST  THE  SECOND  FILE  NAME 
C 

ENTRY  HEABR2 

28.  CALL  SYMBOL  (4. 5, 3-5, *15 , ’FILE  ’,.0,5) 
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P'l.  CAM.  SYMBOL  (A.  ,3.5  , .  IS  ,IFILF.P(  l),.0,ft) 

BO.  RCTifRN 

U.  END 


1.  SUBROUTINE  FTAXES 

2.  LOGICAL  TWOPLT 

3.  COMMON  /FTCOMM/  PLTMIN,PLTSTP,AXS, ISP, TWOPLT 

U.  CALL  FACTOR ( 0 . 3  ) 

CALL  PL0T(1. ,.5,-3) 

CALL  AXISM(0. ,0. ,'DB' ,2,R./AXS,90. .PLTMTN ,AXS*PLTSTP, . IS , .15 ) 
CALL  AXIS(0.  ,0.  ,'FREQ/  ,-U,10.  ,0.  ,0.  ,.05*ISP) 
r 

C  DRAW  SECOND  AXES  IF  TWO  PLOTS,  ELSE  RETURN 

r 

IF  { .NOT. TWOPLT)  RETURN 

C 

CALL  PL0T(.0,»*.5,-3) 

CALL  AX ISM ( .0,.0,'DB' ,2,8./AXS,OG. ,PLTMIN ,AXS*PLTSTP, . 15 , . 15 ) 
CALL  AXIS(.0,.0,'FREQ*  ,-l»,10.  ,0.  ,0.  ,.05»ISP) 

CALI.  H.OT{.0,-U.5,-3) 

RETURN 
END 


SUBROUTINE  CLOSE( NTAPE,NTAP2 ,DISK) 

C 

C  THIS  SUBROUTINE  IS  SYSTEM  DEPENDENT  AND  IS  NEEDED 
C  IF  ANY  SPECIAL  ACTION  HAS  TO  BE  TAKEN  TO  CLOSE  THE 
0  INPUT  UNIT  OR  THE  TEMPORARY  D1..K  IF  USED. 
n 

C  NTAPE  IS  THE  INPUT  UNIT,  NTAP2  IS  THE  TEMPORARY  DISK 
C  AND  THE  LOGICAL  VARIABLE  DISK  IS  'TRUE’  IF  TEMPORARY 
0  DISK  WAS  USED. 

r 

C  THE  FOLLOWING  CODE  SHOULD  BE  CHANGED  AS  REQUIRED 
C 

LOGICAL  DISK 

DIMENSION  NTAPE( 5 ) ,NTAP?{ 5 ) 

C 

C  CLOSE  AND  CLEAN  TEMPORARY  DISK  IF  USED 
IF  (.NOT. DISK)  CALL  CL0SEl(0 ,NTAP2) 

C 

C  CLOSE  THE  INPUT  UNIT 

CALL  CLOSElfl , NTAPE) 


C 


RETURN 

END 
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SUBROUTINE  BREAD(  IT4PE.IFIL  ,NUM,ARF.A  .NSTART) 

C 

C  THIS  SUBROUTINE  READS  NUM  ELEMENTS  FROM  FILE  'IFIL' 

C  STARTING  AT  ELEMENT  NSTART.  THE  ELEMENTS  ARE  TRANSFERRED 

C  INTO  THE  REAL  ARRAY  ’AREA' 

c 

C  THE  FILE  RESIDES  ON  TAPE  OR  DISK  AS  SPECIFIED  BY  ITAPS(l) 

0  INT  THE  PRESENT  SYSTEM,  THE  MEDIA  IS  AS  FOLLOWS 
r» 

C  ITAPE(l)  *  9998  PERMANENT  DISK  FILE 

0  ITAPE(l)  *  9999  FILE  ON  TEMPORARY  DISK 

(’  OTHER  VALUES  FILE  ON  TAPE 

0 

0  IN  THE  CARE  OF  A  TAPE,  ^HF.  TAP*’  LABEL  IS  STORED  IN 

r  ITAPEO-3).  in  all  cases  the  filename  is  stored  IN  IF  I L  ( 1-2 ) 

c 

C  NOTE.  CALLS  TO  9909  CAN  ARISE  WHEN  READING  A  TEMPORARY  FILE 

G  CREATED  BY  TEKPLOT,  OR  IF  THE  USERS  ORIGINAL  DATA  FILE 

C  RESIDED  ON  THE  TEMPORARY  DISK 
C 

C  THIS  SUBROUTINE  IS  SYSTEM  DEPENDENT,  AND  THE  FOLLOWING  CODE 

C  SHOULD  BE  REPLACED  BY  THF  NECESSARY  DATA  READING  CALLS 

C 

DIMENSION  ITAPE(5),IPIL(S),AREA(1^ 

C  WORDS  ITAPE(5)  AND  IFIL(3-5)  USED  ONLY  IN  THIS  SYSTEM 
COMMON  /M1READ/  I ENTER 
IENTER=1 
ITAPE(5)=0 
C 

C  SEPARATE  9999  FROM  OTHER  CASES 
IF  ( ITAPE( 1 ) . EQ . 9999 )  GOTO  1 
CALL  DMREAD(1 .ITAPE , IFIL, NUM, AREA .NSTART, 0 ) 

RETURN 

C 

C  HANDLE  9999  READS 

1  CALL  DMREAD(0, ITAPE, IFIL, NUM, AREA, NSTART, 0) 

RETURN 

END 


SUBROUTINE  BWRITE{ ITAPE, IFIL , NUM, AREA .NSTART) 

C 

C  THIS  SUBROUTINE  IS  THE  INVERSE  OF  BREAD,  AND  WRITES  NUM 

C  ELEMENTS  FROM  A  REAL  ARRAY  'AREA1  ONTO  A  FILE  'IFIL' 

C  STARTING  AT  POSITION  'NSTART*  IN  THF,  FILE. 

C 

C  AS  IN  BREAD,  ITAPE(l)  WOULD  SPECIFIY  WHERE  THE  FILE  RESIDES, 
C  BlfT  IN  TEKPLOT,  BWRITE  IS  ALWAYS  CALLED  WITH  ITAPE(l)=9999, 
C  IE.  FILES  ARE  ONLY  WRITTEN  ONTO  TEMPORARY  DISK. 


THIS  SUBROUTINE  IS  ALSO  SYSTEM  DEPENDENT,  AN!)  THE  FOLLOWING 
CODE  SHOULD  BE  CHANGED  TO  THAT  APPROPRIATE 

DIMENSION  ITAPE( 5 ) ,IFIL< S ) ,AHFA( 1 ) 

COMMON  /M1REAP/  I  ENTER 

IENTER=1 

ITAPF.C  S  )=0 

IX)  THE  TRANSFER 

CALL  DMRITE(0 ,ITAPE,IFIL,NUM,AREA,HSTART,0) 

RETURN 

END 
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